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ABSTRACT 

This paper presents the results of a Spitzer\YiAC 3 — 8 yum photometric search 
for warm dust orbiting 17 nearby, metal-rich white dwarfs, 15 of which apparently 
have hydrogen dominated atmospheres (type DAZ). G166-58, G29-38, and GD 
362 manifest excess emission in their IRAC fluxes and the latter two are known to 
harbor dust grains warm enough to radiate detectable emission at near-infrared 
wavelengths as short as 2 /xm. Their IRAC fluxes display differences compatible 
with a relatively larger amount of cooler dust at GD 362. G 166-58 is presently 
unique in that it appears to exhibit excess flux only at wavelengths longer than 
about 5 /im. Evidence is presented that this mid-infrared emission is most likely 
associated with the white dwarf, indicating that G166-58 bears circumstellar dust 
no warmer than T ~ 400 K. The remaining 14 targets reveal no reliable mid- 
infrared excess, indicating the majority of DAZ stars do not have warm debris 
disks sufficiently opaque to be detected by IRAC. 

Subject headings: circumstellar matter — infrared: stars- minor planets, asteroids- 
planetary systems - stars: abundances — stars: evolution — stars: individual 
(G166-58, G29-38, GD 362, PC 0235+064)— white dwarfs 



1. INTRODUCTION 

The Spitzer Space Telescope opened a new phase space to white dwarf researchers inter- 
ested in the infrared properties of degenerate stars and their environments. The majority of 
white dwarfs are inaccessible from the ground beyond 2.4 /im due to their intri nsic faintnes s 



combined with the ever- increasing sky brightness towards longer wavelengths (lGlasslll999l ). 



This limits any white dwarf science which aims to study matter radiating at T < 1500 K. 
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Prior to the launch of Spitzer, only one previously published, directed mid-infrared study 
of white dwarfs exists ; an Infrared Space Observatory searc h for dust emission around 11 
nearb y white dwarfs, 6 of which have metal-rich photospheres ( IChary. Zuckerman. fc Becklin 



1999|). 



Owing to the superb sensitivity of Spitzer (jWerner et al.ll2004j ). a Cycle 1 IRAC program 
was undertaken to search for warm dust emission associated with cool, hydrogen atmosphere 
white dwarfs with photospheric metals, the DAZ stars. This paper presents a synopsis of 
the IRAC results, including the detection of 5 — 8 /zm flux excess at 0166-58, 3 — 8 /zm data 
on 029-38 and CD 362, and also includes Oemini 3 — 4 /xm spectroscopy of 029-38. 



2. SCIENTIFIC MOTIVATION 



2.1. White Dwarfs Are Metal-Poor 



The origin and abundances of photospheric metals in isolated white dwarfs has been an 
astrophysical curiosity dating back to the era when the first few white dwarfs were finally 
understood to be sublu minous via the cornbination of spectra and parallax; Sirius B, 40 Eri 
B, and van Maane n 2 (Ivan Maanenlll919l : lAdamslll915l . Il914l ). In a half page journal entry, 
van Maanenl (119171 ) noted that his a ccidentally discover ed faint star with large proper motion 
had a spectral type of "about FO" (Ivan Maanenl 119171 ). It was not until forty y ears ('!) later 
that it became clear that van Maanen 2 was metal-poor with respect to the Sun (jWeidemann 
1960). Over the course of the another decade and a half, it became gradually clear that 
white dwarfs in general - both hydrogen atmosphere degenerates void of metallic features 
and helium atmosphere degenerates with photospheric meta l lines - had heavy element 



abundances a few to several orders of m agnitude below solar (IWehrselll975l : IOrenfelllll974 



Shipman 


1972; 


Weener 


1972; 


Bues 


1970) 



evidently sub-solar, heavy element abundances are now referred to as metal-rich. 

The two fundamental classes of metal-poor- yet-rich white dwarfs are essentially the same 
as their metal-free counterparts, which are designated by the main atmospheric constituent; 
either hydrogen or helium. It is easiest to think of white dwarfs as type DA (hydrogen) or non- 
DA (helium). For a wide range of temperatures, white dwarfs with hydrogen atmospheres will 
manifest Balmer lines (type DA), excepting perhaps the very hottest cases where hydrogen 
may be ionized (e.g. PC 1159 stars). White dwarfs with helium atmospheres display a variety 
of spectral behaviors which depend primarily on their effective temperatures; helium II lines 
(type DO) for the hottest, helium I lines (type DB) for a wide temperature range, and no 
lines (type DC) for the coolest. Although there do exist cases of mixed atmospheres, either 
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inferred or directly observed, typically one of the two light gases dominates the composition 
and photosphere; a notable exception is the case of carbon opacity in a helium atmosphere 
(type DQ). There are both subtleties and complexities in the composition and spectral 
characteristics of white dwarfs which will not be discussed here, but for the full guided tour 
of the white dwarf sp ectral zoo, the re ader is referred to both the definition of the curren t 
classification scheme ( Sion et aD Il983l ) and a broad spectral atlas ( Wesemael et al. 1993 ). 
The important points to take away from a given designated spectral type is: 1) the first 
letter after D indicates the dominant gas in the atmosphere; 2) when metals are detected in 
an ultraviolet or optical spectrum of a white dwarf, the letter Z is added to the designation 
(with DC becoming DZ). 



2.2. Metals Are Contaminants 



Any primordial heavy elements in white dwarfs can only be sustained in their pho- 
tospheres for the brief {t ~ 10^ yr) period while the degenerate is still rather hot and 
contracting significantly, and then only to a certain degree. This is achieved through ra- 
diative levitation at T^ff > 20, OOP K for hydrogen atmo s pheres and T^ff > 30, OOP K for 



helium atmospheres (IChayer. Fontaine, fc Wesemaellll995 



Fontaine fc Michaudlll979f ). Be- 



low these temperatures, the cooling degenerate stars develop significant convection zones. 



(Chaver. Fontaine. & Wesemaell 1995: 


J^aaue 


tte et al.lll986: Muchmore 


1984 


Alcock & Illarionov 


1980; 


Vauclair. Vauclair. & Greenstein 


1979; 


Fontaine & Michaud 1979 


: Fontaine & Van Horn 


1976; 


Schatzman 1958). Diffusion timescales for the sinking: of metals are alwavs orders of 



magnitude shorter than the evolutionary (cooling) timescales of white dwarfs. Therefore, 
external sources are responsible for the presence of the metals within cool white dwarf pho- 
tospheres. 

Although there were suspects and many spurious detections reported, all metal-rich 
white dwarfs were historic ally helium atmosph ere degenerates until the confirmation of 
G74-7, the first DAZ star (ILacombe et al.lll983l ). Even so, a decade later all other known 
metal-bearing white dwarfs were still restricted to t he helium-rich variety, despite r eports 
to the contrary; an examination of Tables 1 and 2 in iDupuis. Fontaine, fc Wesemael ([1993') 
with hindsight reveals that only G74-7 remains class ified as DAZ ( IZuckerman et al.l l2003i : 
Wolff. Koester. fc Liebertll2002l : IWesemael et al.lll993l ). This situation began to change a lit- 
tle less than a decade ago wit h a couple of individual discoveries, followed by many more, now 



total l ing around fifty objects (IKoester et al 



19981 : iHolberg. Barstow. fc Green 



1997 



2005 



Zuckerman et al. 



2003 



Zuckerman fc Reid 



Koester. Provencal, fc ShipmanI 119971 ). The reason 
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DAZ stars are latecomers to the metal scene is because helium atmospheres are quite trans- 
parent relative to hydrogen atmospheres; all else being equal, a given calcium abundance 
will produce an equivalent width of order 10^ — 10^ times stronger in a helium as opposed 
to hydrogen atmosphere, making it far easier to infer the presence of photospheric me tals in 
helium- rich degenerates (IZuckerman et al.ll2003l : iDupuis. Fontaine, fc Wesemaell 119931 ) . 

Convection zones can be many orders of magnitude larger in heliur n atmosphere white 

10^ vr JPuDuis et a l"l992': 



Paauette et al. 


19861; 


Muchmorel 


1984; 


Vauclair. Vauclair. &; Greenstein 


1979 



Alcock fc lUarionovlEgSol: iFontaine fc Michaudil979i : 
Fontaine fc Van Horn 19761 ). This relatively long- 



lived photospheric retention allows for the possibility that extant metals in such objects are 
the remnant of a long ago (up to several diffusion times c ales) interstellar cl oud encounter 
(IDupuis et al.l Il993l : iDupuis. Fontaine. &: Wesemaell Il993l : iDupuis et al.l Il992l ) . However, it 
was shown that the Galactic positions a nd space motions of DZ and DBZ stars are no t 
correlated with local interstellar clouds (lAannestad et al.l Il993l : lAannestad fc Sioru Il985l ). 
Additionally with regards to cool helium atmosphere white dwarfs, there is the ever-present 
and still unexplained lack of detectable hydrogen (or more specifically in the case of those 
bearing metal lines, the very high inferred metal to hydrogen abundance ratios), which should 
be expected in quantity, and readily visible in low opacity helium-dominanted photospheres, 
i f accretion from the interstellar me dium has occurred at any time during their cooling 
(IDupuis. Fontaine, fc Wesemaell Il993l ). 

The timecales for the diffusion of metals in hydrogen atmosphere white dwarfs diminish 
quite drastically relative to those in their helium atmosphere counterparts, owing to sig- 
nificantly smaller convection zones for T > 6000 K, and in warmer cases (T > 13000 K) 
where the convection zone is almost negligible, these diffusion times c an be just a few days 
( iKoester fc Wilkenll2006l : IZuckerman et al.ll2003l : iPaquette et al.lll986l ). Unless such objects 
have just emerged from interstellar accretion episodes within these short periods (unlikely), 
they must be cu r rently accreting at rates sufficient to produce the observed abundances. 



Zuckerman et al.l (120031 ) find no correlation between enhanced int erstellar medium d e nsities 
and the current positions of accreting DAZ stars. In contrast, iKoester fc WilkenI (120061 ) 
conclude that, with certain caveats, the warm, partially ionized medium can explain the 
observed accretion rates and abundances in DAZ stars. 



2.3. Minor Planets Are Metal-Rich 



Two decades ago it was proposed that accretion of circumstellar material might be the 
cause of the heavy metal abundances seen in some white dwarfs, so that the origin of con- 



- 5 - 



taminating elements seen is bimodal in nat ure; either interstellar or circumstellar (jSion et al. 
1990l : lAlcock. Fristrom. fc Siegelmarull986l ). The first such model invoked episodic cometary 
impacts from reservoirs which managed t o survive the post-main sequence mass loss phases 



particularly the asymptotic giant bra nch (jPebes fc Sigurdssonll2002l : iParriott fc Alcocklll998 



Alcock. Fristrom. fc Siegelmanlll986l ). This particular model fails to explain 1) the DAZ stars 
with the highest metal abun dances, and 2) the observed distribution of abundances in gen- 
eral (jZuckerman et al.ll2003l ). A more promising model of circumstellar accretion invokes the 
tidal disruption of an asteroid, which goes on to form a r ing of debr is around the white dwarf, 
from which the photospheric heavy elements originate (|Jurall2003l ). Given that a typical so- 
lar system asteroid is around 10^ times more massive than a typical comet, this model can 
explain both the relatively high metal abundan c e and the observed infrared excess seen at 
sever al metal-rich white dwarfs (jJura et al.ll2007l : iFarihi et al.ll2007l : iBecklin et al.ll2005l : lJura 
20031 ) ■ While it may be the case that both mechanisms create contaminated white dwarf 
photospheres - perhaps circumstellar in the case of high metal abundances and interstellar 
for those in the lower range - there are a growing number of me tal-rich white dwarfs which 
are e i ther confirmed or suspected to harbor circurnstellar dust (fjura. Farihi. fc Zuckermaii 
200?! : Ivon Hippel et al-lboO?! : iKihc et alj[2006l : ljuralbood : iReach et aliboOSah . 



Observations in the mid-infrared are most sensitive to both warm and cool orbiting dust 
at metal-bearing white dwarfs, providing a direct test of circumstellar accretion hypotheses. 
Additionally, such a search can constrain the frequency of orbiting material more strongly 
than ground-based near-infrared observations, which are sensitive only to dust not far from 
its sublimation temperature. 



OBSERVATIONS AND DATA 



3.1. IRAC Imaging of DAZ Degenerates 



Table [Ulists the 17 metal-rich wh ite dwarfs observed with IRAC, taken from lZuckerman et al. 



(120031 ) with the exception of GD 362 (iGianninas. Dufour. &: Bergeronll2004j ) . All targets were 
chosen as DAZ stars, but r ecent evidence impli es that two of the white dwarfs (G77-50 and 
GD 362) are helium-rich (IKoester et al.l l2005l : D. Koester 2007, private communication). 
For the purposes of this paper and the statistics which follow, many targets are nominally 
referred to as DAZ, even though future observations may reveal helium-rich atmospheres 
in some. This fact is simply unavoidable as hehum becomes spectroscopically undetectable 
in white dwarf s cooler than Tpff ~ 10, 000 K and its presence can only be inferred by indi- 



rect methods. (IBergeron. Saffer. fc Liebertlll992l ). These stars were chosen among available 
Spitzer Cycle 1 observations for their high metal abundances. 
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Between 2004 Nov ember and 2005 A ugust, observations were executed with the Infrared 
Array Camera (IRAC: lFazio et al.ll2004l ) in all four bandpasses: 3.6, 4.5, 5.7, and 7.9 fim. A 
20-point cycling dither pattern (of medium step size) was used for each target in each band- 
pass, with 30 s frame times at each position, yielding a total exposure time of 600 s at all 
wavelengths. The data were processed with the IRAC calibration pipeline (versions 10 — 12) 
to create a single, fully processed and reduced image (1".20 pixels in all four channels) upon 
which to perform measurements. Aperture photometry was carried out with the standard 
IRAF task apphot, and measured fluxes were corrected for aperture size, but not for color. 
Generally, the flux and signal-to-noise were measured in a 2 — 3 pixel aperture radius, depen- 
dent upon target brightness and neighboring sources, with a 10 — 20 pixel sky annulus. The 
measured fluxes were converted to the standard IRAC a perture using corrections des cribed 
in the most recent version of the IRAC Data Handbook (jSpitzer Science Centerll2006l ). The 
results are listed in Table [H 



3.2. Photometric Errors 

Listed in Table [2] are the total errors in the measured and calibrated flux, together with 
IRAC pipeline versions with which the data were extracted from the archive. The photo- 
metric errors were estimated by taking the per pixel standard deviation in the extracted 
sky level and multiplying by the area of a 2 pixel radius photometric aperture, the small- 
est radius for which there ex ist derived aperture corrections in the IRAC Data Handbook 



( ISpitzer Science Centeiil2006l ). This approach is conservative and does not assume Gaussian, 
random noise for the following reasons. It was found in general that the observations went 
sufficiently deep as to approach or reach the confusion limit, primarily, but not exclusively, 
at the two shorter wavelengths. Additionally, there was sometimes low spatial frequency 
structure seen in the background of the two longer wavelength images, which appeared to be 
from either real diffuse sources such as cirrus, or from imperfections in the IRAC pipeline. 
These factors led to difficulty in the determination of the true sky level in the vicinity of the 
science target, the per pixel noise in the sky, and in the flux measurement itself. Even in 
the smallest, 2 pixel radius, aperture used for photometry, there was occasional possible flux 
contamination from neighboring sources in this relatively large 18.1 square arcsecond area, 
compounded by the undersampled pixels. This possibihty was evident from a few (typical) to 
several (rare) percent increase in photometric flux that was sometimes seen as the aperture 
radius was increased from 2 to 3 pixels, which should not occur for single point sources and 
accurate aperture corrections. All targets were unambiguously detected at all wavelengths. 



There are flve other sources of photometric error in the IRAC camera which were con- 
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sidered ( ISpitzer Science Centerl 120061 ) . First, the absolu te calibration uncertainty in the 
IRAC instrument is reported as 3% (iReach et al.l l2005bl ) . Second, color corrections have 
been ignored and are typically smaller than 1%. Third, there is the pixel phase dependent 
correction, which is reported to be no rn ore than 4% peak to peak or ±2% for a single image 
at 3.6 fj,m (ISpitzer Science Centeiil2006l ). The images analyzed here were produced by a set 
of 20 pseudo-random dithers which essentially annihilates this source of error. Fourth, there 
is the array location dependent correction for Rayleigh- Jeans type stellar sources. This effect 
is the largest source of photometric uncertainty in the IRAC instrument and may reach 10% 
peak to peak or ±5%, therefore a typical error of this type should be 2 — 3%. However, ex- 
periments have shown that for well dithered data, as th at analyzed here, the effect tends to 
average out (as would be expected) and is less than 1% (ICareyll2006l : ISpitzer Science Center 
2006h . 



To assess the reliability of the absolute IRAC calibration, a couple of experiments were 
undertaken (M. Jura 2007, private communication). First, the 3.6 to 4.5 /im flux ratios of 
the target stars with Tes > 7000 K, and no evidence for the presence of warm dust, were 
examined and found to vary typically within 5% of the mean, however a couple of targets 
deviated by as much as 10% - quite significant for such a small number of white dwarfs. 
Second, a similar result was obtained by exa mining 4.5 to 7.9 /xm flux ratios of identically 
warm white dwarfs from lMuUally et al.l (120071 ). In this instance the rel atively lower signal-to- 
noise at 7.9 fim may have been a f actor as eviden c ed by the result of iTremblay fc Bergeron 
(120071 ) who fit a subsample of the iMuUally et al.l (120071 ) white dwarfs with current models 
to within 4% at 4.5 /im, yet only to within 10% at 7.9 fim. White dwarf models themselves 
cannot account for the observed deviations in flux ratios at these warmer temperatures, which 
can only amount to around 1% for extreme values of s urface gravi t y and temperature (P. 
Bergeron 2007, private communication). Furthermore, iHines et al.l (120061 ) report observed 
IRAC fluxes for 33 main sequence stars which deviate from model predic ted photospheric 
values, on average, by 8% at 3.6 /xm, 4% at 4.5 /im, and 6% at 7.9 fim. ISilverstone et al. 



(120061 ) find similar deviations between measured IRAC fluxes and model predictions for 74 
young main sequence stars, amounting to 5% on average at all wavelengths. Based on these 
analyses and findings, it seems appropriate to assign a Icr IRAC calibration uncertainty of 
5%. For bright targets, the total error is dominated by the absolute calibration uncertainty, 
while for the faintest targets the total error is dominated by the uncertainty in the aperture 
photometry. For these reasons, the total errors represented in Table [2] should be considered 
conservative. 
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3.3. NIRI L-grism Spectroscopy of G29-38 

On three nights in 2006 January G29- 38 was observed at G emini Observatory at Mauna 



Kea with the Near- Infrared Imager (NIRI; iHodapp et al.ll2003l ) in spectroscopy mode using 



the L-grism, which covers 3.0 — 4.1 /xm. Spectra were taken at two positions along the 0".75 
slit with 1 s exposures and 60 coadds. Overall, approximately 2 hours of usable science 
frames were gathered. Calibration frames taken each night included spectral flats as well as 
observations of one of the two AOV telluric standards HIP 110578 or HIP 12640, utilizing 
0.2 s exposures and 10 — 20 coadds nodded at two positions along the slit. 

Both the science and telluric standard frames from each night were pairwise subtracted 
at the two nod positions in order to best remove the bright, variable sky at these wavelengths. 
The subtracted frames were fiat fielded and then median combined, with bad pixels and 
cosmic rays fixed manually, creating two spectra of opposite polarity which were extracted 
using standard IRAF tasks. The two extracted spectra from each night were wavelength 
calibrated using aborption lines from the night sky, and averaged. The science spectrum each 
night was divided by the telluric spectrum and multiplied by a blackbody of the appropriate 
temperature. The science spectra from all three nights were shifted and averaged to create 
the final spectrum. The fiux was converted from Fx to F,^, normalized to one, then fiux 
calibrated using IRAC 3.6 fim photometry. The signal-to-noise was estimated by measuring 
the standard deviation along 20 sections of 50 pixels each spanning the entire wavelength 
range, yielding values between 4 and 10 (average 7) over 3.0 — 3.4 /im, and between 8 and 
15 (average of 12) over 3.4 — 4.1 fim. These were calculated assuming continuum or pseudo- 
continuum over all regions. 



ANALYSIS AND RESULTS 



4.1. Spectral Energy Distributions 



The IRAC fluxes for all targets, together with optical and near-infrared data, are plotted 
as spectral energy distributions in Figures [THSl In general, the short wavelength photometry 



was taken from the r aost a ccurate and reliable sources avail able, including 



to: 



McCook fc Sion 



flDENIS Consortium 



Monet et al.ll2003 



20031) and references therein; 2M ASS (ISkrutskie et al. 



20051): and various other sources (jLiebert. Bergeron, fc Holberg 



3 ut no t hmited 
2006h: DENIS 



Zuckerman et al.ll2003l : iBergeron. Leggett. fc Ruizll200ll : iBergeron. Ruiz. &: Leggett 



!rg6 



2005 



19971 ). Because the aim of the survey was to identify significant mid-infrared photometric 
excess due to opaque dust, blackbody fits to the spectral energy distributions of the target 
white dwarfs will suffice to model their expected, essentially Rayleigh- Jeans, behavior at 
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IRAC wavelengths. If the total uncertainty in the IRAC photometry decreases sufficiently, 
white dwarf models should prove more useful in this regard. 

Of the 17 observed stars, three display excess radiation within their IRAC beams (con- 
sidered here to be the solid angle contained within a full width at half maximum Airy disk, 
typically 2".0 — 2". 4 in diameter and 3.1 — 4.5 square arcseconds in area at 3 — 8 /im), with 
a high degree of certainty. The case of G166-58 is discussed in some detail below and shows 
evidence for continuum emission from T ~ 400 K dust. Both G29-38 and GD 362 show 
warm (T ^ 900 K) thermal continuum and strong silicate emission in t heir mid-infrared 
spectra and have th erefore been conffimed to harbor orbiting rings of dust (jJura et al.ll2007l : 
Reach et al.ll2005al ). For the remainder of the targets, the IRAC observations either rule out 
emission from warm, opaque debris or are of sufficiently low signal-to-noise as to preclude 
a definitive conclusion either way. Unfortunately, there was a single target (G21-16) for 
which accurate fluxes could not be extracted due to unavoidable source confusion within 
and around the IRAC beam. An unsuccessful attempt was made to photometrically isolate 
the white dwarf with the IRAF task daophot; based on the signifiant crowding, it is likely 
that the IRAC beam itself is contaminated at all wavelengths. 

From the present work, 3 of 17 metal-rich white dwarfs, or 18% display definite IRAC 
continuum excess consistent with orbiting dust. However, this fraction is only 2 in 15 DAZ 
stars, or 13%. If one counts all DAZ white dwarfs observ ed with IRAC in Cycle 1, the total 
which display IRAC flux excess is 3 of 25 targets or 12 % (Ivon Hippel et al.ll2007l : iKilic et al. 
2006h . 



4.2. GD 362 



The spectral energy distribution of GD 362, together with its IRAC dat a, is displayed 



in Fi gure [51 As expected on the basis of ground-based JHKL'N' photometry (IBecklin et al. 



20051 ) . the 3 — 8 /xm region is dominated by warm thermal emission from orbiting dust, whose 
flux can be reproduced out to A ~ 6 /xm by a T = 900 K blackbody. This approximation 
is physically nonviable as it corresponds to a ring (or sphere) with a single temperature 
and ra dius, yet does not deviate much from a ring model of substantial radial extent lJura 
(120031 ) at these shorter mid-infrared wavelengths. The measured flux at 7.9 /im has not been 
correcte d for color and shows an indication that the broad and strong 10 /im silicate emission 
feature (jJura et al.ll2007l : iFarihi et al.ll2007l ) contributes significantly into that relatively wide 
(AA = 6.5 — 9.5 /im) bandpass. 



In Figure El the 2 — 24 /im photometric data on GD 362 are plotted together with the 
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more realistic model of I Jural (120031 ) . The model invokes a face-on geometrically thin opaque 
dust ring of finite radial extent, with an inner temperature of Tin = 1200 K and an outer 
temperature range Tout = 300 — 600 K. For GD 362, these ter aperatures correspond to a 
ring which extends from ~ 0.1 Rq to Dout ~ 0.3 — 0.7 Rq (IChiang fc GoldreichI 119971 ). 
The data agree reasonably well with the higher temperature curve, excepting the 24 /im 
flux (disregarding the silicate emission-enhanced 7.9 /xm flux). There are two possibilities 
based on this model: 1) the outermost opaque grains have a temperature Tout < 600 K 
or; 2) the flux at 24 /im is affected by another dust emission feature (which would indicate 
forsterite, if present). Photometric and spectroscopic observations of GD 362 utilizing all 
three instruments aboard Spitzer, including a detailed model fit of the orbiting dust, its 
mid-infrared thermal continuum and emission features, dust mass estimates, temperature 
and particle size distribution, dimensions and types of the emitting regions are presented in 



Jura et al.l (120071 ). 



4.3. G29-38 

The spectral energy distribution of G29-38, together with its IRAC data, is displayed 
in Figure [5l While the overall similarity between G29-38 and GD 362 is apparent, there are 
some distinctions. Although their inner dust temperatures are clearly similar, the thermal 
continuum flux of G29-38 over 3 — 6 /xm appears to be falling, while for GD 362 it is rising; 
likely an indication of varying amounts of opaque dust both slightly warmer and cooler than 
the T = 900 K blackbody approximations for the excess at each of these stars. In order to 
fit the slope of the 3 — 6 /xm photometry of G29-38 with a blackbody, a temperature near 
1100 K is necessary, greatly overpredicting the near-infrared flux. Yet as can be seen from 
the figure, the 3.6 /im flux is somewhat underpredicted by 900 K. 



Using the more plausible and physical model of I Jural (120031 ) . Figure [6] plots the thermal 



infrared excess of G29-38 from 2—24 /im from all available Spitzer photometric d ata, togethe r 



with the model curves. These disk models are exactly those applied to G29-38 in I Jural (120031 ). 
now plotted with more accurate data with greater wavelength coverage. For the same inner 
and outer temperatures given above, the extent of the opaque ring is from Di^ ^ 0.1 Rq to 
-Dout ~ 0.4 — 0.9 Rq for G29-38. In this case, unlike GD 362, all the photometric data are 
fitted decently by the model where the outermost grains have a temperature near 600 K, 
with the possible exception of the 3.6 /im flux. 

The apparently excessive flux from G29-38 at the shortest wavelength IRAC channel 
was first noticed when the data were initially retrieved from the Spitzer archive using pipeline 
version 11.0. It was present again one version later, and then triple checked in 2006 November 
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with version 14.0, which was used for all the IRAC data on G29-39 in this work. The 
deviation between the models and the measured flux at 3.6 fim is 0.92 — 0.98 mJy from the 
best fit ring model or the blackbody. Taking photometric error into account, which is entirely 
due to instrument calibration uncertainty, the average deviation is 0.70 mJy; approximately 
2.8(7 of the flux error, or 10% of the excess. 3 — 4 /im spectroscopy was undertaken to 
investigate possible sources of this extra emission, and to better assess if it is real. 

As can be seen from Figure [TJ the L-grism spectrum of G29-38 is essentially featureless, 
with an apparent, slight, upward slope towards 4 /im (a thermal continuum approximated 
by a 900 K blackbody would peak longward of 5 /im in F^). Previous to this investiga- 
tion, the only spectral informatio n in th is wavelength regime came from data presented in 



Tokunaga. Becklin. &: ZuckermanI (jl990l ). where a relatively low signal-to- noise spectrum is 
not inconsistent with emission over 3.2 — 3.7 /im. The spectral flux error in the Figure [7] 
data, translated to mJy via the IRAC 3.6 /im flux, is typically 1.2 mJy over 3.0 — 3.4 /im 
(outside the atmospheric transmission window, an area very sensitive to water vapor and 
prone to large variability), and 0.7 mJy over 3.4 — 4.1 /im (atop the L'-band). 

There exist a plethora of emission features in this region associated with polycyclic aro- 
matic hydrocarbons and their multitudinous c lose relatives: specifically, features between 3.2 
and 3.6 nm seen in the i nterstellar medium (IDraine 



1989; 



2003; 



AUamandola. Tielens. &: Barker 



Geballe et al.lll985l ): circumstellar matter (IMalfait et al 



Geballe et al 



1998; 



Beintema et al.l 



1996 



1985[) ; ultrayiolet -excited (planetary, p roto-planetary, and reflection) nebulae 



fiGeballe et al.lll992l. 



19951: 



19891. 



19851): as well as in comets (IBockelee-Moryan. Brooke. &: Crovisier 



Brooke. Tokunaga. fc Knacke 1991; Baas. Geballe. fc Walthei 19861 ). The well known 



3.3 /im feature, when present, is a lways significantly weaker than i ts other family mem- 
bers at 6.2, 7.7. 8.6. and 11.3 am JPrainel lioosl: iMalfait et al1ll998l : iBeintema et~d1ll996 



AUamandola. Tielens. fc Barked 1 198 91). Because these s tronger features are absent from the 



IRS 5 — 15 /im spectrum of G29-38 (IReach et al.ll2005al ). it is expected no feature should be 
present at 3.3 /im, although the IRS data were published after the the L-grism observations 
were planned. There are additional features around 3.4 /tm (primarily due to methanol, 
ethane, and oth er hydrocarbon species) which typically dominate this region when observed 



towards comets (Mumma et al 



2001 



^ . Crovisier et al.lll997l : lBockelee-Morvan. Brooke, fc Crovisier 

I995I : iBaas. Geballe. fc Waltheil Il986[). but which a re much weaker than the 3.3 /im feature 
in circumstellar environments (IMalfait et al.l Il998l : iBeintema et al.l Il996l ) . Because of these 
observational facts, in addition to the relatively fragile nature of hydrocarbons (some small 
specie s are more volatile than water ice) i r i the vicinity of high density ultrav iolet radiation 
fields ( jJoblin. Boissel. fc de Parsevallll997l : iBrooke. Tokunaga. fc Knackdl 19911 ) , it is unlikely 
that a cometary feature would be seen in the vicinity of the 0.1 — 0.4 Rq circumstellar dust 
ring at G29-38. 



- 12 - 



The failure to detect any possible sources of excess emission in the 3 — 4 /xm region 
at G29-38 leaves a few possibilities for the 2.8 a disagreement between its measured IRAC 
flux and the applied models. The first is that the photometry in the 3.6 /im bandpass is 
impr o perly calibrated, making it inaccurate. The second is that both the ring model of 



Jural (120031 ) and the single temperature blackbody fail to predict the correct flux at this 
wavelength. A T = 1000 K blackbody (Figure [6]) is not inconsistent with the 2 — 6 /xm 
photometry, but underpredicts the three longer wavelength Spitzer data points, requiring 
that dust emission affects those bandpasses. The third possibility is the discrepancy arises 
from photometric variability. It is well-known that G29-38 is a pulsating white dwarf; i?-band 
light curves reveal periods of 615, 268, 243, and 186 s with amplitudes of 0.12, 0.03, 0.03, and 
0.02 mag respectively, which are mirr ored at -ft^-band with corresponding streng ths of 0.02, 
0.02, 0.03, and 0.03 mag respectively (jPatterson et al.lll99ll : iGraham et al.lll990l ). Although 
matching periods were searched for at L, none were found, but with sufficiently large uppe r 
limits which do not exclude variations similar to those seen at K (jPatterson et al.lll99ll ). 
Because the IRAC observations lasted 600 s, these photometric variability timescales cannot 
expl ain the 3.6 /xm flu x being possibly high. Given that the IRAC channel 1 (AA = 3.2 — 3.9 
fim; 
/xm; 



Fazio et a. 



Tokunaga 



2004r) data agree reasonably well with ground-based L-band (AA = 3.2 — 3.8 



20001) data with no color corrections to either, the discrepancy is likely to be 



model disagreement. The 4.5 and 7.9 /im fluxes f or G29-38 in Table P agree well with the 
quoted (uncorrected for color) fluxes reported by iReach et al.l (l2005al ). It is possible that 
the deviation arises from variability which has not yet been seen or reported, but no such 
claim is being made based on the present data. A fourth possibility might be a companion, 
but all substellar models of appropriate age range predict about twice as much excess at 4.5 
fim relative to 3.6 /im, as well as a commensurate rise from 3 to 4 nm which is no t observed 
in the L-grism spectrum (IBaraffe et al.ll2003l : iBurrows. Sudarsky. fc Lunind 120031 ). 



4.4. G166-58 



Figure m displays the spectral energy distribution of G166-58 together with its IRAC flu x 
measurements. The plotted data are opti cal BVRI from [B erger on. Leggett. fc Ruia (120011 ) . 
the average of 2 [/-band values cited in iMcCook fc Sion (20 03il. and near-infrared JHK 



which are the average of values given in lZuckerman et al.l (120031 ) and lBergeron. Leggett. fc Ruiz 

( 2001 ). Th e effective temperature o f the white dwarf is taken from both lLiebert. Bergeron, fc Holberg 
(12005! ) and IZuckerman et al.l (120031 ). where the determinations are both within 0.6% of 7400 
K. G166-58 displays clear excess emission within its IRAC beam beginning at 5 /xm, making 
it unique among white dwarfs conflrmed or suspected to harbor orbiting dust. In all other 



cases, the excess becomes unambiguous by 3 /xm, as in G29-38 and CD 362 (Ivon Hippel et al. 
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20071 : iKilic et al.ll2006l : iBecklin et al.ll2005l : IZuckerman fc Becklinlll987bl l The implied tem- 
perature of the 5 — 8 yum excess can be reproduced by a blackbody of T ^ 400 K, and is 
substantially cooler than the T ^ 900 K temperatures inferred in all other dusty white dwarfs 
with 2 — 3 /xm excess. Yet circumstellar dust at 400 K is still considered warm relative to 
the overwhelming majority of main sequence stars with infrared excess attributed to a debris 



and Dr] 



)h„. > 10 AU 


(Beichman et al. 


2006; 


Sons; et al. 


2005; 


Laureiis et al. 


2002; 


Zuckerman 


Chen & Jura 


2001 


)• 



Before proceeding further with any analysis and interpretation, owing both to the 
uniqueness of the IRAC data on G166-58 and the IRAC field in its vicinity, the nature 
and validity of the excess must be examined. 



4.4- i- IRAC Beam Contamination 

There is another source 5". 3 from G166-58 seen in all four IRAC images, which are dis- 
played in Figure [HI This source is also present in the SDSS Photometric C atalog (release 5) 



where it is designated SDSS J145806. 96-1-293726. 3 and classified as a galaxy (lAdelman-McCarthy et al 



20071 ). This object has optical and mid- infrared colors that are consistent with an extragalac- 
tic source. The IRAC images of G166-58 and the nearby galaxy overlap near 2".5 — 3".0 from 
their image centers. Both sources appear point-like at all four wavelengths. The percentage 
of flux from the nearby galaxy in a 2 pixel radius centered at G166-58, measured relative to 
the total flux of the white dwarf over the same area, is relatively benign at 2 — 3%, with the 
exception of 7.9 fim where it is 12%. This was determined by folding the IRAC point spread 
function of 0166-58 in two along the axis of symmetry (up-down in the Figure [H] images) , 
and subtracting one side from the other. Exploiting this symmetry, it was straightforward to 
remove the contamination in the aperture photometry of G 166-58 and vice versa for similar 
measurements of the nearby galaxy (its flux is 0.18, 0.17, 0.13, and 0.43 mJy at 3.6, 4.5, 5.7, 
and 7.9 fim respectively). The values in Table [2] for G166-58 correspond to the signal after 
removal of this unwanted flux. 

Two additional methods used to obtain the flux of G 166-58 are discussed below, and 
a comparison of the fluxes obtained for the white dwarf by each method is shown in Table 
[31 First, the white dwarf and nearby galaxy were photometrically fitted and spatially de- 
convolved using daophot. Second, photometry was obtained by radial profile analysis of the 
white dwarf. Figure [9] displays overlaying and identical linear contour plots for 7.9 /im images 
of G166-58 both before and after the removal of the nearby galaxy. The parameters from a 
2.0 pixel gaussian radial profile fit (image centroid and full width at half maximum) of the 
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white dwarf remain essentially unchanged by inclusion of the galaxy. Radial profile analysis 
of G166-58 at each IRAC channel yields a measurement for the white dwarf, via comparison 
with the radial profiles and fluxes of G29-38 and GD 362. Examination of G166-58 in a 7.9 
fj^m image where only the galaxy is fitted and subtracted reveals a clear point-like source 
and excess at the location of the white dwarf. In summary, Table |3] illustrates excellent 
agreement among the three photometric methods used to determine the flux of G166-58, 
with the exception of daophot at 3.6 /im (which may be due to undersampled data, which is 
most germaine at this shortest wavelength.) 

Further support that the two IRAC sources are separated photometrically is evidenced 
by their 5.7 to 7.9 /xm colors; the nearby galaxy has a flux ratio of 0.30, while the excess 
detected at G166-58 has a flux ratio of 0.46. Hence, with regards to the adjacent galaxy, 
there should be no doubt of its lack of influence in the white dwarf data presented. 

The known density of background galaxies, along with the presence of the relatively 
bright nearby source on the IRAC chips, prompted an evaluation of the probability that yet 
another, hidden source could be contaminating flux measurements within the beam of G166- 
58 itself, and therefore be responsible for the 5 — 8 /im excess. The white dwarf is located 
at galactic latitude b = +68°, hence any contaminating source would almost certainly be 
extragalactic in nature. To assess the probability of such a line of sight coincidence, the 
following estimations were made. The 7.9 /xm excess emissi on at G 166-58 is ab out 0.07 



mJy, or 14.9 mag. Spitzer IRAC 7.9 /im source counts from iFazio et al.l (l2004bl ) indicate 
approximately 3,000 galaxies per magnitude per square degree at 15*^* magnitude. Taking 
the distribution to be flat over a 0.5 magnitude interval centered at 15.0 mag, this yields 
around 1500 galaxies per square degree of appropriate brightness to reproduce the 7.9 /im 
excess seen at G166-58. Because the white dwarf and its excess display a point-like nature at 
this wavelength (see Figure [H]), any unresolved background source would have to lie within a 
small fraction of the IRAC beam width of G166-58 in the plane of the sky, certainly within 
an area of 2 square arcseconds, conservatively speaking. Therefore, the probability of finding 
a galaxy of 0.07 mJy brightness at 7.9 /im within 2 square arcseconds of G166-58 is about 
1 in 4300. The odds that 1 in 17 target stars is contaminated by such a background galaxy 
should then be around 1 in 250. 



Therefore, the most likely explanation for the excess 5 — 8 /im emission at G166-58 is 
circumstellar dust associated with the white dwarf. 
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4- 4- 2- Circumstellar Dust 



If the 5 — 8 fim IRAC fluxes are the sum of the white dwarf and another source, one 
can immediately rule out a cold companion. A substellar object with an energy distribution 
similar to that implied by the 400 K blackbody fit, would have a radius of roughly 2 x 10^ m, 
more than twice the size of Jupiter. Furthermore, the combined spectral energy distribution 
appears nothing like what might be expected from s uch an orbiting cold degenerate; it 



most notably lacks significant flux e xpected at 4.5 fim (IFarihi. Zuckerman. fc BecklinI 12005 



Burrows. Sudarsky. fc Luninell2003l ). Therefore, the excess emission must be due to warm 
circumstellar material. 

As a first stab at modeling the excess, one might invok e tidal dust rings siniilar t o 
those which have been suc cessful for both G29- 3 8 and GD 362 (jjura et al.ll2007l : |Jurall2003l ). 
Following the formalism of lChiang &: Goldreichl (119971 ) , a fiat opaque ring (or disk) geometry, 
implies a dust grain temperature-radius relation given by 



. gr 



3/4 



(1) 



where Tgr, D, Tes, and R are the temperature of the emitting grains, their distance from 
the star, the stellar effective temperature, and the stellar radius, respectively. If one as- 
sumes that G16 6-58 is a single, carbon-oxygen core white dwarf with log g = 7.97 and 
Tpff = 7390 K Jliebert. Bergeron. fc HolbergI [2005[). then its radius is R = 0.0132 Rq 
( iBergeron. Saumon. fc Wesemaell 119951 ). These stellar parameters yield a distance of Din = 
0.38 Rq to 400 K grains. If one assumes the ring extends to where the grain temperature 
is 200 K, then the ring would extend to D^ut = 0.97 Rq. These inner and outer disk radii 
would be significantly larger than those implied for G29-38, where the ring probably does 
not extend much further than -Dout ~ 0.4 Rq (around 30 stellar radii) for outer grain tem- 
peratures of Tout = 600 K, which seem to fit the data in Figure [6] q uite well. This is also 
true of the implied size of the tidal ring about GD 362, where models (I Jura et al.l 120071 ) also 
yield an outer radius Z^out ~ 0.4 Rq (around 40 stellar radii). 

Regarding the outer edge of the disk an d an appropriate scale for the tidal breakup of 
a rocky body such as an asteroid or comet, iDavidssonI (119991 ) provides a thorough review 
and revision of effective Roche limits. The distance, S, at which a small orbiting body will 
be disrupted by the gravitational field of a large body of radius R, can be expressed, in 
simplified form, as 
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pX 1/3 

6^a{-\ R 



(2) 



where P and p are the densities of the large and small bodies, and a is a constant which 
typically has a value in the range 1 — 2.45, but can be smaller. The coefficient a depends 
on the model, which may include factors such as; composition, heterogeneity, size, shape, 
rotation, orbital characterisitcs, shear and tensile strengths. The classical value of a = 2.45 
is for the case of a rotating, uniform, self-gravitating liquid in a circular orbit, whereas for 
a nonrotating, spherical satellite of solid rock or ice, the value becomes a = 1.26. This last 
approximation is lik ely to be valid for typical asteroids and comets with radii r > 5 km, at 
least to within 50% (jDavidssonlll999l : iBoss. Cameron, fc Beng|l99ll ). Taking p = 1 g cm~^, 
which lies within the ran ge of densities for both asteroid s and comets (past = 1-0 — 3.5 g cm^^. 



0.1 — 1.1 g cm" 



-3. 



Binzel. Hanner. &: Steell l2000l ) . and calculating the average density 
of G 166-58 with the parameters above, an estimate of the Roche limit using a = 1.26 is 
6 ~ 1.2 Rq. Hence, if an optically thick, fiat ring orbits G166-58 at distances corresponding 
to 200 K dust grains, these would lie within a region consistent with tidal disruption of a 
minor planet. However, at grain temperatures of 100 K, the implied distance from G166-58 
would be 2.4 Rq, a region where asteroids or comets should remain intact, even if liquified 



(!Rochelll848l ). 



Integrating the fiux of the 400 K blackbody fit to the IRAC excess yields Ljr = 1.0 x 10^^ 
Lq, while the stellar luminosity given by models for a hydro gen atmosphere white dwarf with 
Tpff = 7400 K and log g = S, yields L = 4.4 x 10""^ Lq (IBergeron. Saumon. fc Wesemael 
19951 ). Together, these determine r = Li^/L = 0.0023, which is logarithmically about 
midway between r ^ 0.0002 fou nd for the main sequence A3 star C Leporis, and r 



0.03 for both G29-3 8 and GD 362 flJura et al.l 120071 : iBeckhn et al.ll2005l : iReach et al.ll2005al : 
Chen fc Jural l200ll ). Therefore, instead of invoking an optically thick, fiat ring, it may 
be more appropriate to suppose an op tically thin shell o f blackbody grains in radiative 
equilibrium, whose distance is given by (IChen fc Jurall200ll ) 



T ^ 



R \ 

2D ^''^ 



(3) 



In this case, 400 K grains would be located near 2.3 Rq, about 6 times further out than 
predicted by the opaque disk model and well beyond the Roche limit for large rocks. An 
advantage of this assumption is that it allows an estimation of the minimum dust mass 
contained in the disk. Because white dwarfs have masses M ~ 1 Mq, yet greatly reduced 
luminosities L ~ 10~^ — 10~^ Lq, radiation pressure on dust grains cannot compete with 
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gravitational attraction. This can be seen by examining the ratio of these two forces, repre- 
sented by the parameter 

B = — m 

^ IQ-nGMcpa ^ ' 

where a is the du st particle radius and is the radiation pressure couphng coefficient 



( lArtymowica Il988l ) . Assuming the case of geometric optics, where the effective grain cross 
section equals its geometric cross section and ^ 1, yields the maximum possible value for 
/3 near a = 0.1 fim. For p = 1 g cm'^, L = 4.3 x 10"^ Lq, M = 0.58 Mq, /J^ax = 0.004 and 
hence gravitational forces dominate over radiation pressure. This simply serves to show that 
sub-micron size dust and gas particles could certainly be present at white dwarfs without 
any danger of being lost to radiation pressure. The minimum dust mass of an optically thin 
disk is approximately 

Mdust ~ — j^D pa (5) 

Although both smaller and larger particles are almost certainly present, the size of the 
thermally emitting dust is on the order of 1 — 10 /im. Taking the density of silicate grains to 
be 2.5 g cm~^ gives a rough lower limit to the mass of the the dust disk of Mjust ~ 2 x 10^^ 
g, a mass of a very large comet. While radiation pressure cannot remove dust at white 
dw arfs, drag forces can. The tim escale for Poynting-Robertson removal of particles is given 



by (IBurns. Lamy. fc Soteiill979l ) 



3 LQpj- 

For white dwarfs in general, the ratio D'^/L for optically thin dust will be 10^ — 10^ times 
smaller than for main sequence stars. Silicates of 1 pm size and p = 2.5 g cm~'^ orbiting 
G166-58 at 2.3 Rq will be removed by the drag force in tpr = 460 yr. The ratio Qpr/a 
increases to a maximum near a = 0.1 pm, implying removal timescales about 10 times 
more rapid. For smaller particles, the radiation coupling efficiency swiftly declines and this 
ratio decreases, lev elling off near a = 0.01 pm at a value similar to that for 1 pm grains 



(lArtymowicg Il988l ) . Hence, in an optically thin disk, all particles up to 1 pm in size should 
be removed by Poynting-Robertson drag within 500 yr. Assuming a present balance between 
accretion and diffusion with the above stellar parameters for G166-58, the rate at which it 
currently gains refractory circumstellar material is M = 2.0 x 10® g s"^ (Koester & Wilkens 
2006; where a factor of 0.01 has been included to reflect the absence of accreted hydrogen 
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and helium). If this accretion occurs over a single diffusion timescale of 10^'^ yr, the total 
accreted mass would 8 x 10^*^ g, equivalent to the mass of a small solar system asteroid. 



4-.4-3. Double Degeneracy 



Based on available data, it appears possible or perhaps likely that G 166-58 is not a 
single white dwarf. Although this target was not listed in Table 2 of IZuckerman et al. 



( 120031 ) for white dwarfs known or suspected to be in binary systems, they determine a low 
surface gravity of log g = 7.58 via a combination of optical through near-infrared pho- 
tometry and parallax. At 7400 K, contemporary white dwarf models predict that such a 
low surface gravity implies a radius R ^ 0.016 -Rfr^, which is about 25% to o large for a 
normal carbon-oxygen core degenerate (IBergeron. Saumon. fc Wesemaellll995l ). An inferred 
oversized radius and corresponding overluminosity can be explained either by a single, low 
mass, helium core degenerate or via binarity involving two components of similar (Am < 1 
mag) brightness - a double degenerate. Actually, the fact that G166-58 appears overlumi- 
nous based on its spectral e nergy distribution and parallax has existed since the analysis of 
Bergeron. Leggett. &: RuizI (120011 ). There, using an essentially identical procedure, a low sur- 
face gravity of log g = 7.66 w as determined, carrying the s ame implication. The near-infrared 
JHK measurements of both IZuckerman et al.l (120031 ) and IBergeron. Leggett. fc Ruia (120011 ) 
agree to within their respective errors. 

Because the implied overluminosity of 0166-58 relies fairly heavily on its trigonometric 
parallax(es), a careful literature search was performed in order to assess all available as- 
trometric data. The value employed by IBergeron. Leggett. fc RuizI ( l200ll ) comes from the 
most recent version of the Yale P arallax Catalog, and is given as tt = 0".0289 ± 0".0041 
( Ivan Altena. Lee, fc Hoffleitlll995l ). The source of the Yale catalog parall ax is the US N aval 
Observatory (H. 0. Harris 2006, private communication), first publis hed in lRoutlyl (119721 ) and 
then updated (and possibly revised) in iHarrington &: DahnI (Il980l ). This parallax is given 
asTT^iative = 0".0277 ± 0".0041 and vTabsoiutc = 0".0298 in their Table 1. iMcOook fc SIot 



( I2OO3I ) quote a value of vr = 0".028 while citing iRoutlyl (119721 ). but this is merely the rela- 
tive (i.e. measured) parallax quoted above. Paying careful attention to detail, one can see 
the Yale and US Naval Observatory values are not identical, but reflect slightly different 
astrometric corrections to obtain the absolute parallaxes, with the Yale galactic model being 
more recent and likely more reliable. Therefore, the only existing parallax measurement for 
0166-58, upon which its overluminosity rests, is vrtrig = 0". 0289 ± 0". 0041 01 d = 34.6l|;3 pc. 



Bergeron. Leggett. &: Ruia (1200 ll ) use this parallax for 0166-58, together with its optical 



and near-infrared spectral energy distribution to deduce T^q = 7310K, M = 0.41 Mq, and 
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My = 12.9 mag. However, fitting model s to th e slope and Balmer line profiles of an optical 
spectrum, iLiebert. Bergeron, fc Holbera (120051 ) determine a nearly identical temperature 
but log g = 7.97, M = 0.58 M©, My = 13.28, and d = 29.1 pc instead. This disparity 
still can be explained by the presence of another white dwarf which contributes the extra 
flux at V, but not by a single low mass white dwarf with an overly large radius, which 
is inconsistent with the Balmer line profile fit. The apparent discrepancy might also be 
explained by supposing that G166-58 is located close r to 29 pc, which is nea r the 30.3 pc 
lower limit implied by the uncertainty in its parallax. IZuckerman et al.l (120031 ) remark that 
radial velocity measurements in 1998 June, 1999 April and July agree within their errors and 
hence there is some weak evidence against radial velocity variability in G166-58. A DA or DC 
white dwarf companion at AV ~ 1 mag would dilute the Bal mer lines and possibly redden the 



optic al spectrum, depending on its effective temperature (IBergeron. Greenstein. &: Liebert 



1990l ). It may be the case that such a compani o n sho uld have already caused notice in the 



data analyzed by lLiebert. Bergeron, fc HolbergI (120051 ) . but nothing was noted. 

Figure H] appears to imply that a single temperature blackbody does not fit the optical 
and near- infrared data perfectly. In fact, there seems to be a slight near-infrared excess 
at JHK, which would become more prominent if a higher temperature model were used 
to fit UBVRI only. Yet the IRAC 3 — 4 /xm data lie very close to the plotted blackbody, 
so perhaps there is another explanation for any apparent mismatch between model and 
published flux, such as data which are not photometric or calibration errors. In any case, 
follow up observations of G166-58 would help to evaluate the possibility that it may be 
a double degenerate. Specifically, a radial velocity study, a careful model analysis of its 
spectral energy distribution with accurate and precise photometry, or another trigonometric 
parallax measurement would all be useful. 



4- 4- 4- Circumbinary Debris 



Identification of a double degenerate suspect follows more or less as it does for main 
sequence stars; the object lies superior to its expected position in a Hertzsprung-Russell (or 
equivalent) diagram. This requires that the distance and effective temperature of the star are 
known or constrained in some fashion. For white dwarfs in the field, the nominal sequence 
is located near and about a radius corresponding to M = 0.6 Mf^, or log q = 8 fo r cool to 



warm white dwarfs (IBergeron. Leggett. fc Ruizll200ll : [Bergeron. Liebert. fc Fullbrightlll995c ; 
Bergeron. Saffer. &: LiebertI 19921 ). A single target which lies above the log g = 8 sequence 
can either be a single white dwarf responsible for the bulk of the luminosity via an overly large 
radius (and a comparably low mass and surface gravity), or a near equal brightness binary. 
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Somewhat oxymoronically, single, low mass (M < 0.45 M©, helium core) white dwarfs are 
understood to be the end products of close binary evolution and more likely than not still 



attached to their stellar cannibal s (IHansen &: Phinneyl Il998l : iMarsh. Dhillon. fc Duckl Il995 



Bergeron. Saffer. &: LiebertI Il992l ). But a single white dwarf with low surface gravity will 
have appropriately th inned Balmer lines due to a reductio n in Stark broadening, and this is 
not seen in G166-58 (ILiebert. Bergeron, fc Holbergll2005l ). Therefore, if G166-58 is binary, 
it must be composed of two relatively normal mass white dwarfs of similar brightness. 



Given the fact that several such double degenerate suspects have turned out to be bona 



fide ( 



2001 



arihi. Becklin. &: Zuckermanl 20051 : IZuckerman et al 



20031 : iBergeron. Leggett. fc Ruiz 



Marsh. Dhillon. fc Ducklll995l ). it is appropriate to consider dust models which conform 



to this distinct possibility in addition to those above for a single white dwarf (see Table H] for 
a summary of possible parameters). The average stellar parameters for G166-58 from the two 
analyses which utilize its trigonometric p arallax, both of which f i nd a similar overluminosity, 



are 



og g = 7.62, and T^g = 7340 K (IZuckerman et al.l l2003l : IBergeron. Leggett. fc Ruiz 
20011 ). The resulting effective radiati ng surface for such a solution is R = 0.0163 Rq 



(iBergeron. Saumon. fc Wesemaellll99,'j ). Applying the opaque, flat ring model and Equa- 
tion (1) with these stellar parameters yields a distance of Din = 0.47 Rq to 400 K grains. 
This is a valid model for a single, low mass white dwarf with the above large radius, but not 
for a binary. In short, it is not possible to fit two similar white dwarfs within such a tight 
disk while avoiding interactions (mass transfer) and at the same time maintaining relatively 
cool, and sufficiently distant 400 K grains. Additionally, there would be gravitational inter- 
actions between the binary and such a disk which would probably preclude a flat geometry. 
Returning to the optically thin case. Equation (3) places 400 K grains at 2.7 Rq. With some 
adjustment of para meters, this scenario allows a mple space for a double degenerate to orbit 
without interaction (IMorales-Rueda et al.l 120051 ). and sufficient distance from the binary in 
order to maintain 400 K dust. 

Observations of G166-58 with the Fine Guidance Sensors aboard the Hubble Space Tele- 
scope show that it is spatially unresolved to approximately 0".008 (E. P. Nelan 2006, private 
communication). This precludes separations wider than 0.3 AU or 60 Rq, yet still permits 
binarity with dust located at one white dwarf component. 



4.5. Metal-Rich Double Degenerates 

In Table [H there are three confirmed or suspected double degenerates: G77-50, EC 
1124—293, G166-58. These white dwarfs are relisted in Table [5] together with their divergent 
spectroscopic and photometric parameters. In light of the potential binarity of G166-58, a 
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brief focus on the other similar systems is appropriate. 



G77-50. This white dwarf has a recent parallax measurement from lSmart et al.l (120031 ) 
of TTtrig = 0".0595 ± 0".0032 over a 6.2 yr baseline, placing it firmly at (i = 16.8lQ;g pc. 
Thus, its derived spectroscopic parameters in Table [5l which consist of a good Tgff deter- 
mination from optical and near-i nfrared photometry, but only a c rude log g estimate from 
a rough fit to a weak Ha feature (IBergeron. Ruiz, fc Leggettl 119971 ) . are not consistent with 
the astrometric distance. As can be seen in the Table, the relatively low spectroscopic 
log g value makes G77-50 appear underluminous; this is atypical as overluminosity is the 
sig npost of binarity. This dis c repari cy is resolved by presuming the Ha profile examined 
by [Bergeron. Ruiz. &: Leggettl (119971) is due to two ve l ocity- shifted hues (indeed the line 
profile shown in Figure 23 of IBergeron. Ruiz. &: Leggettl (119971 ) appears asymmetric) or per- 
haps weak magnetism which also causes Balmer line broadening and mutation in cool white 



dwarfs where the lines are already wea k (IZuckerman et al.ll2003l : IBergeron. Leggett. fc Ruiz 
200ll : IBergeron. Ruiz, fc Leggettl 1 19971 1. Its binarity has been almost certainly confirmed via 
the detection of two Ca lines with disparat e velocities, observed a t two epochs with both 
lines revealing individual velocity varation (IZuckerman et al.ll2003l ). Therefore, G77-50 is 
a double degenerate in which both components are metal-rich, one or both of which have 
detectable hydrogen lines. For equally luminous components, models predict M = 0.93 Mq 
(log g = 8.53), and My = 15.74 mag for Tes = 5200 K hydrogen atmosphere white dwarfs. 
Since both components are polluted with metals, only accretion from circumbinary or inter- 
stell ar material is consist ent with the observations and a circumstellar origin can be ruled 
out. iKoester et al.l (120051 ) lists this system as helium-rich based on the absence of H/3 (D. 



Koester 2007, private communication) 



EC 1124— 293. The parallax reported in IBergeron. Leggett. fc Ruia (120011 ) is an un- 
published trigonometric measurement (M. T. Ruiz 2006, private communication) of TTtrig = 
0".0164 ± 0".0017 which implies a very low mass white dwarf or a binary (or both), if ac- 
curate. The difference in the absolute magnitudes of the spectroscopic [My = 12.39 mag) 
and photometric [My = 11.09 mag) parameter determinations implies an extra source with 
My = 11.48 mag which is obviously much brighter than an equally luminous companion to 
the spectroscopically identified star. Therefore, either the system contains a very low mass 
DZ white dwarf with a proportionally large radius so that it dominates the binary spectral 
energy dist ribution, or the trig; onometric parallax is inaccurate. Given the spectroscopic 
analysis of iKoester et al.l (120011 ). which yields a very normal 0.6 Mq DA white dwarf, it 
is difficult to cur r ently reconcile this system as binary. Radial velocity measurements by 
Zuckerman et al.l (120031 ) in 1998 December and 1999 April agree within the errors and give 
a gravitational redshift corrected velocity of Vr = 1 km s~^ for the H/3 line (no errors are 
given, but are likely no greater than a few km s~^), while similar measurements reported by 
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Pauli et al.l (120061 . l2003l ) yield f,. = — 3 ± 4 km s ^ (no epoch given). 



Zuckerman et al.l (120031 ) fi nd narrow to damning evidence o f binarity at some known or 



suspe cted double degenerates (IBergeron. Leggett. fc RuizI l200ll : iBergeron. Ruiz, fc Leggett 
19971 ) such as: broad H/? cores (e.g. G141-2, Case 2), presumably from two velocity-shifted 
yet unresolved cores; two separate H/3 cores (e.g. G271-115, G77-50); and a single, variable 
H/5 core (e.g. LHS 1549). Two of these have been confirmed as binary by other methods: 
G141-2 has been spatially resolved with the Hubble Fine Guidance Sensors (E. P. Nelan 
2006, private com munication) and radi al velocity monitoring of LHS 1549 has determined 
its orbital period (INelemans et al.ll2005l ). 



4.6. PG 0235+064 



The IRAC photometry of this target was problematic due to a nearby M dwarf common 
proper motion companion which is reported here for the first time. The companion, PG 
0235+064B, is separated from the white dwarf primary by 7". 4 at 344° in the IRAC 3.6 /^m 
image (epoch 2005.6), shown in Figure [TOl Examining archival images from 2MASS, and 
the DSS reveals separations and position angles between the A and B components which 
remain essentially constant between 1950 and 2001. Blinking the 1950 and 1990 DSS frames 
clearly shows the pair moving togeth er over 40 yr . The USNO-Bl.O catalog has yU = 0".18 



yr ^ at ^ = 153° for the white dwarf (IMonet et al.ll2003l ). which would have caused the pair 



to separate by almost 10" over the last 55 yr if the secondary were background. Hence the 
pair is bound. 

PG 0235+064B has reliable 2MASS photometry consistent with an early M dwarf which 
has very likely contaminated some previously published photometry and spectroscopy of the 
white dwarf, even as far as the blue optical region, ca using it to appear too cool (red) 
(IZuckerman et al.ll2003l : iGreen. Schmidt, fc LiebertI 119861). Figure [J shows Tpff = 15,000 K 



(DA3.4) provides a better fit to the white dwarf data (IBergeron. Wesemael. &: Beauchamp 
19951 ) than previously published, lower effective temperatur es corresponding to DA4.4 and 
DAS (IHomeier et al.l 1 19981 : iGreen. Schmidt, fc LiebertI Il986l ) . The shorter wavelength pho- 



tometric data in Figure [T were selected as to be minimally, or not at all contaminated by th e 
cool companio n: B from iGreen. Schmidt, fc LiebertI (119861 ) ; J from ISkrutskie et al.l (120061 ) ; 
and HK froni iKilic et al.l (120061 ). Assuming log g = 8, very close to that determined by 



Homeier et al. (jl998l ) from a spectrum likely to be contar ninated by the M dwarf, the white 
dwar f would lie at c? = 70 pc. Using the Hubble GSC2.2 (ISpace Telescope Science Institute 
200ll ) blue and red magnitudes, one can estimate B ^ 17 mag and B — K ^ 6 for the 
red dwarf. This corresponds to a spectral type near M3 and agrees reasonably well with 
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the expected absolute magnitud e of Mk = 6.7 mag at the estimated white dwarf distance 



flKirkpatrick fc McCarthvlll994l ) 



For all of the IRAC images, the task daophot was used in an attempt to remove the 
light of the M dwarf, but it was found that this generally oversubtracted its flux in the 
region of the white dwarf, and combined with the 5 — 6 mag difference in brightness between 
components, proved unreliable. Instead, the symmetry of the point spread function was 
exploited to self-subtract the flux on the opposite side of the M dwarf at the location of 
the white dwarf, introducing an additional error component equal to the square root of the 
percent flux removed. 



5. DISCUSSION AND CONCLUSIONS 



Although a survey of only 17 stars does not allow robust statistics, it is clear that 
the majority of DAZ white dwarfs do not harbor warm dust of sufficient emitting surface 
area to be detected with IRAC. If most or all of these stars do host circumstellar material, 
the fractional luminosities must be relatively low compared to currently known dusty white 
dwarfs. This could result from a modest amount of dust (as in the zodiacal cloud), large 
particle sizes, or cooler material further from the star. 

Another possibility is that any warm dust produced within the Roche limit of a white 
dwarf is swiftly destroyed through mutual collisions, not unlike ice in a blender, as it orbits 
with Keplerian velocities near 0.003c. In optically thin disks, particles with orbital period p 
will collide on a timescale given by tcou = p/t- The tidal rings at G29-38 and GD 362 have 
been modeled to extend from approximately 0.2 — 0.4 Rq, where a typical orbital period 
is only p = 0.6 hr and the resulting collision timescale for r ~ 0.01 is tcou = 2.5 dy. If 
a sizeable fraction of dust produced in a tidal disruption event is initially optically thin, 
then both collisions and Poynting-Robertson drag will compete to quickly annihilate this 
material. The ratio of these two timescales for dust particles orbiting a distance D from a 
star of mass M can be written as 



7 = ^ = 693 ^^"" (7) 

where M and L are in solar units, D is in AU, p is in g cm~^ , and a is in microns. This 
fraction reaches a minimum for 0.1 /im grains with 1 g cm~^ at the inner disk edge, and 
yields 7min > 10 for all possible white dwarf disk parameters. Table [6] gives minimum values 
for 7 at the inner edges of the disks at G29-38, GD 362, and G166-58. Therefore collisions 
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will erode dust grains faster than they can be removed by angular momentum loss. This 
is also true in the case of optically thick disks where: 1) the bulk of material is shielded 
from starlight, and hence the Poynting-Rob ertson effect is diminished, and 2) the collision 
timescale is less than half the orbital period (lEspositdll993l ). Hence, for a wide range of disk 
densities, it is plausible that mutual collisions within an evolving dust ring at a typical white 
dwarf will result in the relatively rapid self-annihilation of the micron size grains required to 
radiate efficiently at 3 — 30 /im. 

Following the tidal disruption of an asteroid, if one models the dust produced as a 
coUisional cascad e, the expected particle size distribution behaves classically as n{a) oc a~^^^ 
(jPohnanyil Il969l ) . For dust at main sequence stars, this distribution is reshaped on short 
timescales as sub-micron size grains are removed by radiation pressure which, as shown 
above, does not apply in the case of white dwarfs. In the absence of radiation pressure, the 
average particle in will have a size a = 5/3 a^i^. For practical purposes, at white dwarfs one 
can assume a^[^ = 0.01 /xm, where particles are already inefficient absorbers and emitters 
of infrared radiation, and anything smaller approaches the size of gas molecules and atoms. 
With such a distribution of extant dust, 99.7% of the particles will have sizes a < 0.1 fim, 
leaving a paltry fraction of larger particles which could effectively support infrared emission 
from the disk. Clearly, such small particles and gas might be present at most or all white 
dwarfs which show signs of circumstellar accretion such as the DAZ stars. 



On the other hand, the persistence of warm dust disks at several white dwarfs (I Jura. Farihi. fc Zuckerm; 



20071 ) must be explained despite the fact that in some or most cases where it is produced, 
it may also be efficiently destroyed. One possibility is that the disk density (which could 
contain gas) becomes sufficiently high as to damp out collisions in the disk and also make it 
optically thick, thus somewhat protecting it from self-erosion and drag forces simultaneously. 
The evolution of such a dense, fluid-like ring is then dominated by viscous forces (differen- 
tial rotation an d random motions) which cause it to spread, losing energy in the process 
( lEspositdll993l ). The maximum lifetime of such a ring occurs at minimum viscosity 



2 2 

p w p 



where w is the radial extent of the ring and a is the surface mass density (lEspositdll993l ). If 
the mass of a large solar system asteroid, 10^^ g, were spread into a tidal ring of negligible 
height {h < 10 m), a radial extent 0.2 — 0.4 Rq, consisting of micron size particles orbiting 
a typical white dwarf, the resulting volume mass density (0.55 g cm~^) would be sufficiently 
high that the mean free path of particles is on the same order as their size. This could 
effectively damp out collisions, thus minimizing viscosity, and with a resulting surface mass 
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density of cr ~ 550 g cm"^, permit a potential disk lifetime - in the absence of competing 
forces - longer than the Gyr white dwarf cooling timescales. However, for sustained accretion 
rates as low as M = 10^*^ g s~^, a 10^^ g disk would become fully consumed within several 
Myr. 



Large rocks and colder material orbiting at D > 100 Rq will be unaffected by any of 
the aforementioned processes, and such a reservoir of material is strictly necessary to supply 
some fraction of DAZ white dwarfs with photospheric metals, regardless of circumstellar dust 
production (collisional versus tidal) and evolution (persistence versus destruction). 

The overall number of white dwarfs with remnant planetesimal belts may be rather 
high based on a growing number of detections. If one takes 12% (§4.1) as the fraction of 
DAZ stars with circumstellar dust as o bserved by Spitzer to d ate, 20% as the fraction of 
DAZ stars among cool DA white dwarfs ( Zuckerman et al. 2003[). and 80 % as the number of 
cool DA stars among all white dwarfs in the field (lEisenstein et al.ll2006l ). then a lower limit 
to the number of white dwarfs with asteroid-type belts is at least 2%. This fraction could 
be as high as 20% if the majority of metal-rich white dwarfs harbor circumstellar matter, 
which raises important questions about the implied frequency of planete simal belts around 
main-sequence stars and the current detection rate (see the Appendix of |Jurall2006l ). 



Owing to their low luminosities, white dwarfs which may have been polluted by heavy el- 
ements in winds or transferred material from substellar companions (lDebesll2006l : iDobbie et al. 



20051 : IZuckerman et al. 
tions 



2005 



20031: 



Sion fc Starrfiddlll984j) are easily identified with IRAC observa- 



200?!: Hansen. Kulkarni. fc Wiktorowiczll2006l : lFarihi. Becklin. fc ZuckermanI 



■Mullallv et al.ll 

Farihi. Zuckerman. fc BecklinI 2005 ) down to T dwarf temperatures. There is no evi- 



dence of such companions in the data presented here, ruling out all but the coldest brown 
dwarfs, active planets and moons as close orbiting, companion-like polluters (Farihi et al. 
2008, in preparation). 
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Table 1. DAZ^ White Dwarf Targets 



WD 


Name 


Teff (K) 


V (mag) 


[Ca/H] 


References 


0032-175 


G266-135 


9240 


14.94 


-10.20 


1,2 


0235+064 


PG 


15000 


15.5 


-9.03 


1,3 


0322-019 


G77-50 


5220 


16.12 


-11.36 


1,4 


0846+346 


GD 96 


7370 


15.71 


-9.41 


1,5 


1102-183 


EC 


8060 


15.99 


-10.43 


1,5 


1124-293 


EC 


9680 


15.02 


-8.53 


14 6 


1204-136 


EC 


11500 


15.67 


-7.72 


1,7 


1208+576 


G197-47 


5880 


15.78 


-10.96 


1,8 


1344+106 


G63-54 


7110 


15.12 


-11.13 


1,8 


1407+425 


PG 


10010 


15.03 


-9.87 


1,9 


1455+298 


G166-58 


7390 


15.60 


-9.31 


1,8,9 


1632+177 


PG 


10100 


13.05 


-10.75 


1,9 


1633+433 


G180-63 


6690 


14.84 


-8.63 


1,8,9 


1729+371 


GD 362 


10500 


16.23 


-5.1 


7,10 


1826-045 


G21-16 


9480 


14.58 


-8.83 


1,8 


1858+393 


G205-52 


9470 


15.63 


-7.84 


1,5 


2326+049 


G29-38 


11600 


13.04 


-6.93 


1,5 



^G77-50 and GD 362 are helium-rich flKoester et all l2005l : D. 
Koester 2007, private communication 



20031: (6) iKoester et al 



Bergeron. Leggett. fc Ruiz 



2001 



References. — flllZuckerman et al.ll2003l: f2MMermilliodlll986l: (3) 
This work; (A) ISergeron. Ruiz, fc Leggettlll997l: (5) iMcCook fc Sion 



(7) ISalim fc Gouldl l2003l: f8) 



200ll : (9) iLiebert. Bergeron. &: Holberg 



20051 : (10) D. Koester 2007, private communication) 
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Table 2. IRAC Fluxes for White Dwarf Targets 



WD 




-^4.5Atm (/^Jy) 




-p7.9Mm (yWJy) 


Pipeline 




Q«n -1- 1 Q 
ouu HZ lo 


91 zL + 11 


IzLD + 17 

l^U ZL 1 / 


71 + 1 R 

1 1 ZIZ lo 


1 1 n 




119 + 7 




'^fi + 1 7 
ou zn 1 ( 


9Q + 1 R 

zn 1 o 


1 9 4 


0322-019 


543 ± 27 


381 + 19 


253 + 18 


149 + 18 


12.4 


0846+346 


310 + 16 


198 + 10 


107+17 


78 + 20 


10.5 


1102-183 


242 ± 12 


137 + 8 


95 + 18 


63 + 22 


11.0 


1124-293 


350 ± 18 


200 ± 10 


142 ± 16 


77+15 


12.4 


1204-136 


164 + 8 


101 + 6 


64 + 13 


25 + 14 


12.4 


1208+576 


597 ± 30 


367 ± 19 


211 + 18 


129 ± 15 


11.0 


1344+106 


558 ± 28 


372 + 19 


266 + 20 


137 + 20 


12.4 


1407+425 


292 ± 15 


159 + 8 


113 + 15 


73 + 16 


12.4 


1455+298 


357+ 18 


222 + 11 


189 + 17 


155 + 18 


12.4 


1632+177 


1683 + 84 


1049 + 53 


679 + 37 


403 + 26 


11.4 


1633+433 


912 + 46 


623 + 31 


389 + 24 


232 + 19 


11.4 


1729+371 


380 ± 19 


395 ± 20 


425 ± 26 


644 ± 34 


12.4 


1826-045 


714 ± 72 


414 ± 47 


257 ± 33 


164 ± 39 


12.4 


1858+393 


201 ± 10 


116 + 6 


73 + 17 


54+17 


10.5 


2326+049 


8350 ± 420 


8810 ± 440 


8370 ± 420 


8370 ± 420 


14.0 



Note. — Error calculations, including both photometric measurements and instru- 
mental uncertainties are described in §3.2. 
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Table 3. Flux Determinations for G 166-58 



Method 












357 ± 18 


222 ± 11 


189 ± 17 


155 ± 18 




410 ± 23 


236 ± 11 


172 ± 13 


164 ± 20 


3^ 


350 ± 35 


225 ± 23 


170 ± 17 


150 ± 15 



Note. — Sec §4.4.1 for a detailed description of the methods used to 
account for and eliminate any and all flux from the nearby galaxy. 

^Aperture photometry. 

^Point spread function fitting (daophot) photometry. 
'^Radial profile analysis. 
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Table 4. Possible Disk Parameters for G166-58 



Disk Type 


log g 


-Reff (-Rq) 


-Din, thick {Rq) 


-Din,thm {R&) 


Circumstellar 


7.97 


0.0132 


0.38 


2.3 


Circumbinary 


7.62 


0.0163 


0.47^ 


2.7 



Note. — In the case of a wide binary where dust orbits a single, 
normal mass, white dwarf component, the dust is circumstellar. A 
single, lo w mass white dwarf of large radius is ruled out by spec- 
troscopy (ILiebert. Bergeron, fc Holberg||2005l ). 



close binary surrounding by a flat optically thick disk is not 
physically possible (§4.4.4). 



Table 5. Confirmed or Suspected Double Degenerate DAZ Systems 



Star 




Spectroscopic*^ 






Photometric*^ 




References 


Teff (K) 


log g 


d (pc) 


Teff (K) 


log g 


d (pc) 


G77-50 


5220 


7.5 


23.3 


5200 


8.04 


16.8 


1,2,3 


EC 1124-293 


9550 


8.04 


33.6 


9440 


7.10 


61.1 


4,5 


G166-58 


7390 


7.97 


29.1 


7340 


7.62 


35.5 


4,6,7 



^Parameters derived from model fits to spectroscopic continuum flux and Balmer line profiles. 
^Parameters derived from model fits to photometric fluxes and trigonometric parallax. 



Reference s. — (1) Bergeron. Ruiz, fc Leggett 19971; (2) Leggett. Ruiz, fc Bergeron 19981: (3) 



2003: (A) Bergeron. Leggett. fc Ruij 



Liebert. Bergeron, fc Holberg 



2005 



2001 



(5) 



Koester et al. 



2001 



(6) 



Zuckerman et al. 



Smart et al. 



20031 : (7) 
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Table 6. Minimum 7 Values at Inner Disk Edges 



Star 




M(M0) 


log(L/Le) 


An (Re) 


Tmin 


G29-38 


0.030 


0.69 


-2.64 


0.14 


19 


G166-58^ 


0.0023 


0.58 


-3.36 


0.38 


12 


G166-58= 


0.0023 


1.2 


-3.18 


2.7 


30 


GD 362 


0.030 


0.75 


-2.90 


0.12 


34 



Note. — The minimum value of 7 is achieved, realistically, 
at p = 1 g cm^ and a — 0.1 /im. 

'^Single white dwarf with an optically thick circumstellar disk. 
^Double white dwarf with an optically thin circumbinary disk. 
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Fig. 1.— Spectral energy distribution of G266-135, PG 0235+064, G77-50, and GD 96 (see 
§4.1). For PG 0235+064, the light of its nearby M dwarf companion has been removed (see 
§4.6 and Figure fTOl) . and a superior fit to the data is achieved for an effective temperature 
of 15,000 K, contrary to previous, significantly lower estimates. 
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Fig. 2.— Spectral energy distribution of EC 1102-183, EC 1124-293, EC 1204-136, and 
G197-47. 
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10.00 



0.10 
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10.00 : 



T = 7100 K 



0.01 
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..I. 
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Fig. 3.— Spectral energy distribution of G63-54, PG 1407+425, PG 1632+177, and G180-63. 



-43- 



10.00 F 



1.00 r 



0.10 



0.01 




G166-58 



T = 7400 K 

bb 

T = 400 K 

bb 



1.0 10.0 
A (/xm) 



100.0 



Fig. 4. — Spectral energy distribution of G166-58. 
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Fig. 5.— Spectral energy distribution of GD 362, G21-16, G205-52, and G29-38. The IRAC 
photometry of G21-16 is contaminated at all wavelengths. 
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Fig. 6. — Infrared excess flux from GD 362 and G29-38 after subtraction of the ex- 
pected wtite_dwBx£_pl^^ at these wavelengths. The data are K- 



band (ISkruts 



(IReach et al. 



de et al. 2006: Becklin et a. 



20051 



I), IRAC 3 - 8 /xm (this work ), IRS 16 /xm 



2005ar). an d MIPS 24 /xm flJura et al.ll2007l : iReach et al.ll2005ar ). Plotted are 



models from lJural (120031 ) for disks with an inner temperature of 1200 K and outer temper- 
atures of 600 K (dotted lines) and 300 K (dashed lines) respectively. Also shown are 900 K 
and 1000 K blackbody curves (dashed-dotted hues). 
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Fig. 7. — L-grism spectrum of G29-38 taken with NIRI, normalized, then flux caUbratcd 
using the IRAC 3.6 fim photometry. There appears to be a gentle slope towards 4 fim, but 
the spectrum is otherwise featureless. The resolution is A/AA ^ 500 and the data are neither 
binned nor smoothed. 
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3.6 yum 
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■/ 

WD 
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WD 
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h 
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Fig. 8. — IRAC images of G166-58 in all four channels, each spanning a 54" x 54" field 
of view. The top of each image corresponds to position angle 125.7°, increasing counter- 
clockwise. The adjacent extragalactic source (discussed in §4.4.1) is located 5". 3 from the 
white dwar f at position angle 87.6°. Th e proper motion of G166-58 (/x = 0."64 yr~^ at 



9 = 165.2°; iBakos. Sahu. fc NemethI [20021 ). is nearly perpendicular to the direction toward 
the background source. 
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Fig. 9. — Linear contours for 7.9 /im images of G166-58, drawn from 1.45 to 2.80 MJy sr~^ 
in steps of 0.15 MJyr sr~^. The dotted lines show the contours for the central part of the 
image shown in the lower right panel of Figure [HI which includes the nearby galaxy. The 
solid lines show identical contours for a similar image where the galaxy has been fitted and 
removed as described in §4.4.1. The solid contour near offset (1,2) is residual signal after 
subtraction of the galaxy. 



-49- 



PG 0235+064 
3.6 jLtm 



■ 




Fig. 10. — IRAC 3.6/xm image of PG 0235+064 and its M dwarf companion, separated 
by 7" A at 344°. The top of the image corresponds to position angle 253.3°, increasing 
counterclockwise. The white dwarf is quite faint at m^sum = 16.0 mag, while the red dwarf 
dominates with ms.ejum = 10.7 mag. 



